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The aim of this paper is to adapt to certain parabolic partial differential equations an 
expansion method of solution developed by S. Faedo ^ for hyperbolic equations. In order 
to make possible a moderately compact presentation, the equations treated are not the 
most general to which the method is applicable, but are the simplest nontrivial relatives of 
the heat equation. Similarly, the boundary values and initial conditions are not the most 
general, but are assumed to be in a canonical form to which others, if sufficiently smooth, 
can be reduced. The method of solution not only shows the existence of a solution, but 
describes a definite procedure for approximating it. Some remarks are made on the possi- 
bility of estimating the error. 



1. Statement of Problem 

In the region S, 0<x<t, 0<^<oo, we seek a 
solution of the equation 



L{u)=u^^—Ut—gu=f 
subject to the initial conditions 



(1) 



(2) 



The functions / and g are continuous in S, together 
with their first two partial derivatives. As a com- 
patibility condition, we require that /(0,0) =/(t,0) =0. 
It is convenient to assume that ^>0. This is not 
at all essential, since a substitution of the form 
u=e^^v transforms (1) into 



Vxx—^t—(g+^)v=fe ^'. 



(10 



The new coefficient ^+X can be made positive in 
as large a part of S as desired, and this is sufficient 
for our purposes. 

This method can be applied to more general equa- 
tions in which Ut has a variable positive coefficient 
and the Ux term is present. However, the equation 
selected is of sufficient generality to illustrate the 
method and results and yet not make the computa- 
tions unpleasantly complicated. In case g = constant, 
the method reduces to the method of expanding in 
eigenf unctions. 



2. The Moment Method ^ 
We seek an approximation to u in the form ' 
Un{x,t)=^Cn,kii) sin kXy 

k 



(3) 
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2 means 2 and 



1 1 means I • 



subject to 

f{L(Un)-J}^mjxdx=0, i=l,2, . . ., n, (4) 
where 



Equation (4) amounts to 



(5) 



k 



where 



f>k,j(t)=- fg{xf) sin jx sin kxdx. 



(J^j(f)=- y/(^jO sin jxdx. 

TT 



The functions Cn,j{t) are thus determined by a 
system of n simultaneous first-order hnear differential 
equations. These equations are already in canonical 
form, that is, the derivatives are solved for in terms 
of the other quantities, and so they have a unique 
solution satisfying (5). Thus the functions Un(x,t) 
are determined in unique fashion. 

We intend to show that the sequence Un(x,t) 
converges to a solution of (1). The method proceeds 
by showing that the Un and certain of their deriva- 
tives are uniformly equicontinuous. Kecall that the 
functions / of a family are uniformly equicontinu- 
ous in a region, provided that e is given greater 
than 0, there is a 5>0 such that \f{p)—J(g)\<i€ 
for any / of the family and any two points p and q 
of the region whose distance does not exceed 8. 
Now, Ascoli's theorem asserts that out of a uni- 
formly bounded and equicontinuous sequence of 
functions a uniformly continuous subsequence can 
be selected. Therefore, by selecting subsequences 
of subsequences a sufficient number of times, we 
shall arrive at a subsequence of the Un which con- 
verges uniformly, together with sufficiently many 
of its derivatives to allow us to make the existence 
proof. This latter is effected by using a Green's 
function and showing that the -u^, for n->oo , come 
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successively closer to satisfying an integral equation 
corresponding to (1), and thus that the limit function 
does satisfy such an equation. 

To prove equicontinuity of our sequences, we make 
use of a criterion of Tonelli ^ as follows : A class of 
continuous functions v(xf) in the region St'-^^xK-k, 
0<t<T, is uniformly equicontinuous in this region 
if (a) for each x and 7, v(x,t) and v(xj) are absolutely 
continuous functions of t and x, respectively, (b) the 
partial derivatives Vx and Vt are summable in St, 
(c) there exists a constant c such that J^vldx<Cc and 
Sv1dx<Cc for all members of the family. Further- 
more any accumulation function of a class satisfying 
(a), (b), and (c) also satisfies (a) and (b). 

Now the functions Un and those partial derivatives 
in which we are interested are easily seen to be 
continuous and to satisfy (a) and (b), and the 
burden of the proof of equicontinuity falls on obtain- 
ing inequalities of the form in (c) . This is done by 
making use of devices similar to those frequently 
employed to show uniqueness, and to recall them, 
we shall first prove a uniqueness theorem. 

3. Uniqueness 

Suppose that two solutions of the system (1) (2) 
exist, all derivatives appearing in (1) being con- 
tinuous in S. Their difference v will satisfy (1), (2) 
with/ replaced by zero. Thus L{v)=0 and so 



fL{v)vdx = f{ Vxx—Vt — gv}vdx = 0. 



(7) 



Now 



and 



• 1 d 



fvtvdx = f^ -77 v^dx=^ -T- fv^dx, 



-li. 



dt 



dt 



fVxxVdx=VxV 



— fvldx= — fv\dx. 



If we put K{t) = fvHx, (7) becomes 

hK'{t)=-f{vl+gv')dx<0, 

since ^>0. Since by (2), i^(0)=0, we see that 
K{t) < . But K is by its nature nonnegative and 
K=Oj which implies that v^O and that the two 
solutions are identical. 

Theorem 0. There is at most one function, con- 
tinuous together with all derivatives appearing in (1) 
which satisfies the system (1), (2). 

4. Certain Bounds 

This section is devoted to obtaining the inequali- 
ties required to carry out the program outlined in 
section 3. 

4.1. Lemma. Let K{x) he a nonnegative junction 
satisfying J or x>0 the inequality K\x)<a^K{x), 
where a>0. Then for a;>0, ■^K{^<^K(fi) + {al2)x. 

7 Tonelli, L'estremo assoluto degli integral! doppi, Annali della Scuola Normale 
di Pisa (1933). 



Suppose that for some positive Xq, V^W^ 
Vi^(0) + («/2)xo. Then there exists Xj su ch that 
VF(^= V^W + W2):ri, but V^(^> V^(0) + W2)x 
for Xi<Cx<Xq. In particular, i^(x)>0 in Xi<Cx<ixQ,y 
and so in that interval 



K'{x) 



<( 



^K{xy 

Integrating from Xi to Xq, we have 



^K{xo)-^K{x,)<(a/2) (xo-x,), 



or 



^iK(xo) < ^K(xO + (a/2) {xo-xO = ^K{^) + (a/2)xo, 

a contradiction. The lemma is proved. 

Corollary. Let K{x) be a nonnegative function 
satisfying for x>0 the inequality 



K\x)<a{x)'ylK{x), 



(8) 



where a (x) is nonnegative and nondecr easing. Then 
for a; > 0, 



^JKix)<^[K(fij + ^xa(x). 



(9) 



Consider any Xo> . For x <Xo, K^ (x ) <a {xo)-yjK(x) . 
By the preceeding lemma, 'ylK(x)<-yjK{0) + ^a(;Xo)x. 
Setting x=Xo and noting that Xq is arbitrary, we get 
the stated result. 

4.2. The bound ^ for fuldx. Let each of the 
equations (4) be multiplied by Cn,j(t) and the results 
summed. There comes out the equation 



f {Uxx—Ut—gu—f]udx=0. 



(10) 



Since 



Suxxudx= — fuldx and futudx=\ -j- fu^dx 



lA 



dt 



(10) may be written 



\ -r fu^dx= — f{ul-\-gu^)dx—J*fudx<ffu\dx. 

(11) 

Set J'u^{x,t)dx=K{x). From (11) and Schwarz's 
inequality, 



W{^)<f\Mdx<-ylfjHx-4fuHx. (12) 



Set 



ma.xslff{x,r)dx=:^F{t). 



(13) 



Then (12) gives K'(t) <2F(t)-y/K(t). Fr om th e corol- 
lary of 4.1 we deduce that -yjKjt) <'yjK{0)+tF{t). 
Now K(0) = fu\x,0)dx=0, and so -yjKit) <tF{t). 

8 In section 4 we shall frequently omit the subscript n from Un. Whenever 
u appears in section 4, it will mean Un. 
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Theorem 1. 



ful(x,t)dx<t'F'{t), 



(14) 



where F is given by (13). 

4.3. The bound Jor jTultdx. Let the equations (4) 
be differentiated with respect to t, multipHed by 
Cn,j{t), and added. There results from this 



f[uxxt—utt—gut—gtu—ft]utdx=0. 



(15) 



Again, 



SuxxtUtdx= — fultdx and futtUtdt=\^ J'ufdx, 



dt 



so putting K(t) = J'ufdx, we get 

^K'(t) = -f(u% +gu^)dx-f(g,u+f,)u4x (16) 

<fi\gM+\ft\)\ut\dx 



We define 

and 

Then 



<^f(\gM + \ft\rdx-Mi^^ 



max \gt(x,t)\-=Gi(t)j 



max ^/fJKx,r)dx= Flit). 



T<t 



(17) 
(18) 



f{\gMMJt\ydx<f{Gl{tW+2G,{t)Wt\+n)dx 



< Gl{t) fuHx-\-2G,{t)^JfuHxff4x+ fjfdx 

< Gmt'F\t)+2GS)tF{t)F,{t) +F?{t) =H\t), 

this last equation being a definiti on of H{t) . 
From (16) we have K\t) <2H{t)^/K{t), and by 4.1, 



's[K(i)<'y/K{0)+tH{t). 



(19) 



Now 



K{0)=fu^(x,0) = f(J^C:,,{0)smkxydx 

=~i:CUO)y=- JlifAxfi) sin kxdxf 

<fj\x,0)dx 

by BesseFs inequahty. Thus Kifi) < F\0) , and from 
(19) we conclude that V^(0 </^(0)+^i?(0. 



Theorem 2. 

ful,(x,t)dx<(nO)+tH{t)y=A\t), 
where F is given by (13), 

H\t) = G\{t)eF\t)+2G,{t)tF{t)F,{t) +Fl{t), 
and Gi and Fi are given by (17) and (18). 



(20) 
(21) 



4.4. The bound of ful^dx. Equation 10 may be 
written 

S{ —ul—'UUt—gu^—fu}dx=0, 
or 

fuldx= — f{gu^+uut+fu]dx<f(\uut\ + \fu\)dx 



< -yl fu^dx ' fuldx-\- -^ fu^dx ' Sj^dx 



< 'yJfF'{t)A\t) + 'yJt'F\t)F\t) 

==tF(t)(A(t)+F{t)). 

Theorem 3. 

fuUx,t)dx<tF(t) (A{t) +F{t))=B(t), (22) 

where A is given by (20). 

4.5 The bound oj Sultt(x,t)dx. We differentiate 
(4) twice with respect to t, multiply by 0'^j-(t) and 
sum, to get 

f{u^xtt—Uttt—gntt—2gtUt—gttU—ftt]utt=0. (23) 

As before, we write 

J '^xxtt^ttdx= J u^ttdXy 

futuUttdx=\j-^ fu'ttdx, 

and (23) becomes, on setting K(i) = J^UudXj 

^K\t) = — f[u%,+gul}dx-~f[2gtUt+gttU+ftt]uttdx 
< f { \2g{iii\ + \gtA^ \ht\ } \ntt\dx. 

It is clear that by the use of Schwarz's inequality we 
can obtain an inequality of the form 



K'{t)<mms), 



(24) 



where J(i) is an increasing and continuous function 
expressible in terms of G'l, F, Fi, and the similar 
functions G2 and i^2 corresponding to gajti. 

We estimate K(^) as follows. Differentiate (4) 
with respect to i and put ^=0, to obtain 

futt{xfi) smjxdx=S{Uxxt{xfi)—gUt{x,0) 

—ft{x,0) } sin jxdx. (25) 

for j = 1,2, . . ., n. Now i^i;(a:,0) is a linear combi- 
nation of sin jx, j=l, 2, . . . 7^, and by (25) its 
first n Fourier coefficients are the same as those of 
U:cxt{x,0)—gUt{xfi)—ft{xfi). Thus by BesseFs in- 
equality, 

fuh (xfi)dx < f { Uxxt{xfi) -gu^{xf)) -Mxfl) } Hx, (26) 

If in (4) one puts i=0, there results 

— S'^tixfi) sm jxdx= fj{xfi) sin jxdx. 
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1 



Two integrations by parts gives 

—fu^xt(xfi)Qijijx dx=ff^x(x,0)mnjx dx, 

and since Uxxti^fi) is a linear combination of sin jx, 
j=l, 2, . . ., n, we see as above, 



fulAx,0)dx<ff.,{xfi)dx. 



(27) 



By using Schwarz's inequality, (27), and the bounds 
previously obtained, the right hand number of (25) 
can be dominated by a continuous function depend- 
ing only on / and g and their derivatives. That is, 
K(0) is bounded uniformly in ti, and from (24) it 
follows that the same is true of K{t), 

Theorem 4. There exists a continuous function C{t), 
determined byf,g and their derivatives ^ such that 



ful,,{x,t)dx<C{t). 



(28) 



4.6. The hound Jor ful^^dx. From eq (15) we 
have on replacing fuxxt^tdx by — fultdx, 

fultdx=^f{utt+gut+gtu+jt]uttdx. 

Because of the bounds previously obtained, it is 
clear that we have the theorem : 

Theorem 5. There exists a continuous Junction D{t), 
expressible in terms of f and g and their derivatives 
such that 



ful,t(x,t)dx<D(t), 



(29) 



4.7. The bound for Sy^txxdx. If equations (4) be 
multiplied by j'^ sin jx and summed, we have 

fulJx= f {ut+gu+f]u^Jx 

< -yj f {Ut + gU-\-fYdx-yl fulJx, 

Now from the previous bounds we can find a con- 
tinuous function Eit) which depends only onf,g and 
their derivatives, such that 



f[u^+gu+f}Hx<E{t), 



(30) 



Thus we have the theorem : 

Theorem 6. There exists a continuous function E{t) 
independent of n such that 



fulAx,t)dx<E{t). 



(31) 



4.8. An integral inequality . Consider the Green's 
function G(x,^) relative to the equation y^^-\-\y=0 
=y(0) =y(T). Specifically 



G(x,^)= 






-1) 

X) 






and its uniformly convergent Fourier expansion is 

TT 1 J 

For brevity, set ' 

then 



Ut+gu+f=^(u), 
i^x—^{u)=L(u)—f. 



(33) 

(34) 



Multiply (34) by G(x,^) and integrate with respect 
to X to obtain 



-um - fHu) G(x,Odx= f { L{u) -f] G{x,&dx, 



Now 

f{L(u)-f}G(x,^)dx 



(35) 



= Z) ^^ f{L(u)-f}^mjxdx, 

TT 1 J 

___£, sin^ {i(ti)-/}sinix&, 

TT n-\-l J 

because of (4). Now L(u) —f=Uxx—^{u) , and 

n 

Uxx= — ^Cn j(t)f sin^a:; hence 
1 

J^Uxx smjxdx=0 
for j>/i, giving 

f{Uu)-f}G(x,^)dx=-- S ^^ f^(u) sin jxdx. 

By Schwarz's inequality, 

\f{L{u)-J}Gix,^)dx\ 



Now 



and 



<- a/S i ' Jj: (fHu) sin jxdxy. (36) 



CO 1 1 

tf{j'-3n' 



X) {f^i"^) sin jxdxf 

00 ^ 

< S iJ^H"^) sin jxdxf<- f^\u)dx. 

1 ^ 

By (30), f^\u)dx<E{t), and so from (36) we can 
conclude 



\f{L{u)-J}Gix,Odx\< 



2 /X 



i^E{t). 



\'2 



Referring back to (35) we have a theorem: 
Theorem 7. 



\-Unm-fHUn)G(x,Odx\< 



V 



2E{t) 



(37) 
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5 Existence Theorem 

Recall now the discussion of eqiiicontinuity and 
the criterion of Tonelli given at the end of section 2. 
Because of the inequalities contained in theorems 1 
through 6, we see that the functions Un, Unx, and Unt 
are uniformly equicontinuous. Thus by Ascoli's 
theorem, we can pick out, by a suitable diagonal 
process, a subsequence {Vm} of {Un} such that 

Vm{x,t)-^u{x,t), 
Vmx{x,t)-^vix,t), 
Vmt(x,t)-^w(x,t), 

the approach being uniform in each case. Because 
of this uniformity, one sees that v=du/dx and 
w=du/di. The derivatives Uxx, y^xt, "^tt exist almost 
everywhere and are summable. We consider the 
expression 

-Vm(^.t) - fHVrn)G(x,^)dx. (38) 

Since Vm=Un^y theorem (7) tells us that the expression 
(38) tends to zero as m— >oo. Also ^{vm)-^<l>(u) 
uniformly, and G is bounded, so we have 

lmi{ -Vmm-fHVm)G(x,^)d^} 

= -ui^,t) - fHu)Gix,Odx=0. 

The integral J'^{u)0{x,^)dx has a second deriva- 
tive with respect to x, namely, —^{u) = —Ut—gu—j. 
Hence m„ exists, is continuous, and satisfies 



-«x.a,«)+«'(^)=0, 



(39) 



which is identical with (1). Since u obviously 
satisfies (2), it is a solution of the system (1), (2) 
in the region St. But T is arbitrary, and a solution u 
of (1), (2) exists in all S. 

We have shown that a subsequence of Un converges 
to the solution u, together with first order derivatives. 
Actually the entire sequence Un so converges. If 
this were not the case, there would be a subse- 
quence Un- such that, for instance, lim Un,9^u at 

some point x,t. From the sequence Unj we could pick 
another subsequence in the manner of this section, 
such that this subsequence has for limit a solution v 
of the problem. By the uniqueness theorem, u^v^ 
and so lim Unj{x,t)=u{x,t), a contradiction. Thus 
we have 

Theorem 8. The sequence Un{x,t) has for uniform 
limit in any St the unique Junction u{x,t) satisfying 
(1) and (2). Furthermore, Unx~^Ux and Unc^Ut 
uniformly in St. 

6. Convergence in Mean of u^kx 

The solution u{x,t) of the system (1), (2) has a 
Fourier series 

u{Xjt)=^\j{t) sinja:, 
1 



which may be differentiated twice with respect to x, 
and the resulting series converges in the mean to Uxx* 
Set 



Vn=^h{t) mijx. (40) 

1 

J\CnAt)-h{t)) and 
(41) 



If (4) are multiplied by 
summed, we get 

f{L(Un)—f}-^ {Un — Vn)dx=0 



Since L{u)=J, we may write (41) as 

f{L{Un — Vn)—L{u — Vn)}^ (Un—Vn)dx=0. (42) 

Recalling that 

L{u) =u,j,—Ut—gu—f+J=Ux:,—^(u) +/, 
we write (42) in the form 



(52 \2 

yj^('^«-^n)J dx 



= f[—HUn—Vn)+J—L{u—Vn)}^2('^n—Vn)dx, 

and Schwarz's inequality gives 



'Y/|5pK-?^«)j dx 



<^!f[-^{nn-Vn)+j-L{u-Vn)Ydx, (43) 

'Now,— ^{Un — Vn)+j= — {Un— Vn)t—g{Un— Vn)-^^\m\- 

formly, and L(;v—Vn) tends to zero in the mean; thus 
the right hand member of (43) tends to zero and so 



f{^2(nn-Vn)\ 



dx-^0. 



From this and an appropriate triangular inequality 
it follows that 






Theorem 9. The second derivatives Unxx of the 
approximating functions tend in the mean to Uxx- 

It is too much to expect that Unxx-^^^xx uniformly, 
or even pointwise. For example, Unxx{0,t)=0, while 
Uxx{0,t)==f(0,t), and this latter does not need to 
vanish identically. 

7. Estimating the Error 

From theorem 9 it is easy to prove that L{Un)-^L{u) 
in the mean. Set L{Un)=fn. Then 

l.m. 

L(Un — u) =L{Un)—L{u) ^f^—f-^O. 



267518—53- 
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Using the method of 4.2, we prove that 

f{Un-Uydx<t'Fl\t), 

where 



Flit) =max^Jf{MT)-Ux,T) Ydx. 

T<t 



(44) 



Theorem 10. The mean square error satisfies the 
inequality 

f{Un-uydx<t'Fl\t), (45) 

where Fl{t) is given by (44), and F^(t)-^0 as n-^co , 

Theorem 10 unfortunately gives information only 
about the mean square of the error instead of about 
the error at a point, which would be desirable. This 
latter will be mentioned below. Apart from this, 
however, theorem 10 gives a rather curious sort of 
information about the m.s.e. It does not give an 
a priori estimate of this error; that is, we cannot say 
that after n steps the error will have a certain bound, 
because we do not know how rapidly Fn tends to 
zero. On the other hand, we do seem to know more 
than the mere fact that the m.s.e. tends to zero. 
We know that if we take n sufficiently large, the 
m.s.e. will be small and we shall know it to be small. 
What we do not know is when we shall know it! 

In section 5 it was mentioned that Uxt i^ ^ sum- 
mable function. Thus from L(u)=f, it follows that 
Uxxt is summable. From the equation L{Un—u) 
— /w— /, we obtain, using the methods of sections 4.3 
and 4.4, the inequality analogous to (22), 

f{{Un-u)xYdx<tF:{t){At{t)+Ft{t)) = Gt'^\ 

wnere F^ and Al are given by (13) and (20), except 
that jn—j is used insted of / in their definitions. 
Now we are unable to say much about the factor A^ 
since it contains terms of the form S{fnt—ftVdx] 
that is, depends on Uxxt and Unxxt- However, it 
appears quite likely in many cases that At will 
remain bounded as n-^<^ . In any case the inequality 

f{(Un'-u),Ydx<G:{t) 



is a valid one, and in a given computation Gl{t) can 
actually be computed. Setting Ln{x,t)=Un—Uj we 
have from Schwarz's inequality 



\^n{x,t)\= {^^Kx{s,t) ds<^xj'{Anx{s,t) ^ds (46) 
<^Jxf{Anx(s,t)yds<^'^^6Jf): 

In an analogous fashion we obtain by integrating 
from TT instead of from zero. 



An{x,t) <^/{T-x)G:{t), 
From (46) and (47) it follows that 



(47) 



An{x,t)<yl^x(ir-x)G:(t). 
Theorem 11. 

\Un{x,t)—u{x,t)\<'J-x{Tr — x)G*(t), 



(48) 



where 



6:(t)=tF:{t){A:{t)+F:it)], 



T<t 

A:(t)-F:{o)-tH:(t), 

{Hmr=G,{t)t'{F:(t)V+2G,(t)tF:{t)+Fl{t) 

+{n(0}^ 



Fl(t) = mEiX^f{U(x,T)-Ux,r)rdx, 

r<t 

For any stage n of the solution, the function G^(t) 
can be explicitly calculated. If the A^ do remain 
bounded, then G^-^0, and we will get a useful 
estimate of the error. 

Los Angeles, April 24, 1952. 
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